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SUMMARY 


The motions of four fully rolled-up vortices representing the vortex 
system trailing behind cruciform wings are studied by theoretical and 
visual-flow methods. The analysis applies throughout the Mach number 
range . 

Equations are developed for the three-dimensional paths traced by 
the vortices behind a cruciform wing banked 45°, and calculations are 
made of the distance behind the wing at which the upper two vortices 
pass through the lower two. It is found that this "leapfrog” distance 
depends upon the lift coefficient, aspect ratio, and span loading of the 
cruciform wing, and that for low-aspect-ratio cruciform wings leapfrog- 
ging may occur within two chord lengths of the trailing edges. 

The various types of vortex motion to be expected throughout the 
angle -of -attack range are considered in some detail, and the interaction 
of the two vortex sheets shed from the cruciform wing is taken into 
account. Results of some water-tank studies are also presented and 
compared with the theory. 


INTRODUCTION 


The downwash behind plane wings has been studied theoretically by a 
number of authors and considerable attention has been given to the 
rolling up of the trailing vortex sheet. The analysis of Kaden 
(reference l) predicted the distance behind the wing at which the sheet 
may be considered to be fully rolled up into two trailing vortices, and 
this work was later used (in reference 2) to demonstrate the usefulness 
of the single horse shoe -vortex approximation for the calculation of the 
downwash behind wings of low aspect ratio. Since at the present time 
cruciform configurations are largely confined to wings of low aspect 
ratio, the rolling up of the trailing vortex sheets is again of major 
concern, and the behavior of the fully rolled-up vortices is again of 
considerable interest . 
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While in the case of the plane wing the vortex sheet became rolled 
up into a vortex pair which simply moved downward at a uniform speed at 
great distances behind the wing, the analogous problem for cruciform 
wings is necessarily more complicated. Instead of the two rolled-up 
vortices, there are now presumably four (one from each wing panel) and 
their induced effects upon one another are such as to produce quite 
intricate paths of motion. The downwash field, of course, may therefore 
become extremely involved. This report is concerned with a study of the 
motions of the four rolled-up vortices and their effect on the downwash 
behind cruciform wings. 
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SYMBOLS 


aspect ratio 
span of one wing (2s 0 ) 
root chord 

lift coefficient of cruciform wing 

/ pi 

lift coefficient of plane wing ( — 

VqS 

section lift coefficient for plane wing 
distance hehind wing trailing edge 

distance hehind wing trailing edge where four rolled-up vortices 
are collinear (leapfrog distance) 

constant of the motion related to distance between centers of 
gravity of two vortex systems (y x + y 2 ) 

parameter which depends on initial vortex positions 





lift of cruciform wing 
lift of plane wing 
free -stream Mach number 


free -stream dynamic pressure — pU 5 ” ^ 

V2 !■ 


parameter used for periodic motion 


f "(4-G) 

kG 
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Q 2 parameter used for aperiodic motion (-Q 2 ) 

R 2 

h(G + b) 

S area of one wing 

s local semi span of one wing 

s Q maximum semi span of one wing 

s ! semispan of rolled-up vortices from one wing 

t time 

U free-stream velocity 

v velocity of vortex center in y direction 

w velocity of vortex center in z direction 

x,y, z right-hand orthogonal coordinates with x axis in stream 
direction 

y Q , z Q initial values of y x and z ± 

z i } 2. vertical distance between vortices 1 and 2 

(See sketch on page 9.) 

a angle of attack of cruciform wing axis 

a ! angle of attack of one wing 

(3 angle of sideslip of one wing 

7 ratio of vortex strengths 

T circulation (positive counterclockwise) 

p mass density 

Subscripts 1, 3, and 4 refer to vortex number. (See sketch on page 5.) 
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GENERAL ANALYSIS 


The phenomenon of the rolling up of the vortex sheet behind plane 
wings was discussed as early as 1907 by Lanchester (reference 3) and has 
since been the subject of a large number of investigations. As a result 
of such investigations, it has been found (reference 2) that, whereas 
the vortex sheet behind wings of high aspect ratio may be considered to 
remain flat at the tail location, a good approximation to the downwash 
several chord lengths behind wings of low aspect ratio is often obtained 
by considering the sheet to be fully rolled up into two trailing vortices. 
In the case of the cruciform wing, as in the case of the plane wing, the 
rate of rolling up of the vortex sheets will depend upon the aspect ratio 
of each wing. It should therefore be possible to obtain a reasonably 
good approximation to the downwash field behind low-aspect-ratio cruci- 
form wings by considering the vortex sheets to be fully rolled upj that 
is, by replacing the two vortex sheets by four discrete vortices. It is 
with this simplified model of the physical problem that the analysis of 
this report is concerned at the outset. Subsequently, it will be seen 
that the theory can be modified to provide a more accurate representation 
of the actual flow field. 

With low-aspect-ratio wings in source-free flow the linearized dif- 
ferential equation for the perturbation velocity potential cp 

(l-M 2 )9 + 9+9=0 (1) 

x 1 xx yy zz 

can "be satisfactorily approximated at all Mach numbers "by the two- 
dimensional equation 

<Pyy + ^ZZ = 0 ( 2 ) 

since (l-M 2 ) cp^ is much smaller than Cpyy and cp zz if j/i-M 2 A is 
small. (See references 4 and 5-) Therefore, the problem of the motions 
of the four rolled-up vortices can first be treated by lateral strip 
theory as a two-dimensional problem in planes perpendicular to the 
flight direction. For the cruciform wing at 45° angle of bank, where a 
vertical plane of symmetry is present, the problem is thus reduced to 
that of the motion of two pairs of vortices with a common axis. This 
problem was first attacked by Love (reference 6) and later by Hicks 
(reference 7)> both in connection with the analogous three-dimensional 
problem of the motion of vortex rings. Unfortunately, neither of these 
authors worked completely the problem of interest here, the former being 
concerned only with the relative paths of the vortex pairs about each 
other, whereas the latter was concerned chiefly with the actual problem 
of the vortex rings. 

Since this paper is concerned with the three-dimensional paths of 
the four rolled-up vortices behind the wing, the paths in transverse 
planes will first be determined, and will then be related to the elapsed 
time or the distance the wing has traveled in the flight direction. 

The analysis will be confined to the case of 45° bank. 
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Determination of Vortex Paths 


Relative motions in transverse plane .- The motion of any system of 
two-dimensional vortex filaments is determined by the induced velocity of 
each vortex due to all the other vortices of the system. It is shown in 
reference 8 that the induced velocity of any particular filament can be 
obtained from the function 

W = -p; zr m r n log (Sn-Sn) ( 3 ) 


where m<n and = y m + iz m . That is, if one writes 

W = 0 + i\|r (4) 

then \|f is a constant of the motion and is analogous to the stream 
function in giving the components of the velocity of the particular 
vortex. Thus, the path of each vortex can be determined from equa- 
tions ( 3 ) and (*0, provided that the motion of one vortex determines the 
motions of all the other vortices of the system. (Such is the case, for 
instance, where all the vortices but one are images of the one.) In the 
present problem, since only two of the four vortices are images, equa- 
tions ( 3 ) and (4) are not sufficient to describe completely the motion of 
the system. However, these equations will be useful in determining the 
paths of two of the vortices relative to the other two. It will be seen 
later that this information is of considerable interest. 


For the cruciform wing banked k^° , the simplified vortex system as 
viewed in the y,z plane consists of two pairs of two-dimensional 
vortices having a common plane of symmetry as shown in the sketch. 

It is apparent from the symmetry of 
the configuration that 



Further, if each wing of the cruciform 
can be considered as a plane wing in 
sideslip (as was demonstrated in 
reference 9 for the calculation of 
load distribution) it follows that 

r i = -r a ; r 2 = -r 4 

since the vortex system far behind a 
plane lifting wing must consist of two 
equal (and opposite) vortices, 
although their positions will depend 
upon the angle of sideslip while 
their strength will depend upon the 



plane of 
symmetry 


y 
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lift. On this basis, it can be concluded that 

r i = r 2 = -r 3 = -r 4 = r 

Actually, since the vortex sheets trailing behind the cruciform wing 
must influence one another during the rolling-up process itself, the 
assumption that the wings act independently is not strictly applicable 
to considerations of the wake as in the present problem. However, this 
complication will be deferred to a later section and the analysis will 
first be developed on the basis of four vortices of equal strength. 
Hence, from this point on, the subscripts will be omitted in reference 
to the vortex strengths which will be referred to simply as plus or 
minus T . 


When the directions of rotation of the four vortices are taken into 
consideration, equation (3) can be written for the present problem as: 


W = 



log (Ci-C 


2) 


- log (C 2 -C 3 ) + log (C3-C4) - 


log (Ci-£ 4 ) - log (Cx-C 3 ) - log (C 2 -C 4 ) 


(5) 


Noting further that C 4 and C3 are simply the_negative complex conju- 
gates of C 2 and C 2 > that is C 4 = C 3 = , equation (5) can be 

expressed as 


W = 


ir 


ir* 

kn‘ 


log (£x-C 2 ) - log (C 2 + C 2 ) + log ( Cx-Ia) - 

log (Ci + Ci) - log (Cx + C 2 ) - l°g (C 2 + Ci) 

(Ci-y(Cx-y I 


log 


L(C 2 +C 2 )(Ci+Cx)(Ci+C 2 )(C 2 +Ci) j 


(6) 


Due to the symmetry (see sketch on page 5), the vortex system may 
also be considered as two groups of vortices with equal and opposite total 
circulation (one group on either side of the plane of symmetry). For 
such a pair of groups, the centers of gravity 1 of the two groups must 
remain a fixed distance apart. (See reference 10.) In other words, 
the center of gravity of each group must move parallel to the plane of 
symmetry. Since here the strengths are equal in magnitude, this may be 
stated as 


y c .g. 


zyiri _ r(yx+y 2 ) 

Efi 2r 


yi+y 2 f , , 

— — — = = constant 


(7) 


1 The center of gravity of a group of vortices is defined as the center 
of gravity of a similar field of point masses, the mass of each being 
proportional to the strength of the corresponding vortex. Negative 
masses correspond to negative circulations. 
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With this information, equation (6) can he expressed in terms of the 
coordinates y and z and simplified to 


W 


= Hi log / (yi-y s ) 2 + (zi-z &) 


^ lby iys [f e + (z 1 -z fi )2] 


} 


Now, since W is purely imaginary, equation (4) is simply 


r 2 

* ' i? 106 


f (yx 
I 4v 


~ya ) + ( z i -z 2 )‘ 


^yiyaf^+Czi-zs) 2 ] 


(8) 


(9) 


At this point, a moving coordinate system is introduced such that one 
pair of vortices always lies on the y axis . This is done by a change 
of variable 


Z l,2 = Z 1 ~ Z 2 

Now, noting that y 2 = f - y 1 , equation (9) becomes 



( 2yi -f) a + Zl>s 2 - 

%!( f -y a )(f s +z 1 , a 2 ) _ 


(10) 


(11) 


Setting ^ of equation (11) equal to a constant then yields an equation 
for the paths of vortices 1 and 2 relative to each other# Thus, 


or 


(2yi-f) g t zi, a 2 

y 1 (f-y 1 )( fS + z lf3 s ) 


= constant 



determines the relative paths of vortices 1 and 2 (and also 3 and once 
the value of the constant G is determined. 


In order to evaluate G of equation (12), one must determine the 
initial positions of the four vortices# It must be realized, of course, 
that initially (i.e., at the trailing edges of the cruciform wings) there 
are not four '•discrete vortices but a cruciform system of vortex sheets. 
The immediate problem, then, is to replace these sheets by four vortices 
properly located in the £(or y,z) plane# For the present, it will 
suffice to observe that the initial locations of the four vortices must 
lie on the straight lines formed by the trailing edges of the wings. 
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Thus, for k 5° bank, 

z i ,2 = yi + y 2 = f when y x = y Q ( 13 ) 

where the subscript o refers to the initial value of y 1 . The actual 
determination of yQ will be deferred to a later section of this report 
in which the physical problem will be considered in more detail. 

With the boundary condition expressed in equation ( 13 ) , equa- 
tion ( 12 ) can be written in terms of the initial value of y-,^ : 



It will be shown later that y 0 , and hence G, is a function of angle 
of attack for the cruciform-wing configuration. 

While equation (l4) does not completely describe the motion of the 
four vortices, it yields some interesting information regarding the types 
of motion to be expected. For instance, if the relative path described 
by equation (l4) is a closed curve it can be concluded that the motion 
is periodic. The relative paths of the vortices are plotted in 
figure 1 for various values of G, and it is found that the motion is 
periodic for G less than and aperiodic for G greater than 1)-. 

This is in agreement with the findings of Love (reference 6 ) who 
showed that the condition for periodic motion is that 

r < 3 + 2 J2~ ( 15 ) 

y P 

where r = — at z, . = 0. That is, r is the ratio of the lateral 

y x ’ 

displacements from the plane of symmetry when the four vortices are 
collinear (when one pair is passing through the other). This type of 
vortex motion was first discussed by Helmholtz in connection with 
vortex rings having the same axis and circulations in the same 
direction (reference 11 ): 

"The foremost widens and travels more slowly, the pursuer 
shrinks and travels faster, till finally, if th<?ir veloc- 
ities are not too different, it overtakes the first and 
penetrates it. Then the same game goes on in the opposite 
order, so that the rings pass through each other alternately." 
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Complete motions in transverse plane .- Since the relative paths of 
the vortices on either side of the plane of symmetry are given by equa- 
tion (l^), establishing the actual path of any one of the four vortices 
will completely determine the motion of the system in the y, z plane. 
Therefore, the actual path of vortex 1 will now be calculated. 


The determination of any of the individual vortex paths requires 
the solution of the following equation: 


z n 




d y n = 


'v. 


n 


f n 


d yn 


(16) 


where v n and w n are the velocity 
components of the nth vortex in the 
y and z directions, respectively. 
If equation (l 6) is evaluated for 
one vortex, say vortex 1, then the 
path of vortex 2 can he obtained 
from equations (7), (10), and (l4), 
and the paths of vortices 3 and k 
are found by symmetry. It will be 
seen that this is considerably 
easier than solving four equations 
of the type given by equation (l6). 
By adding the contributions of all 
the other vortices to the velocity 
components of vortex 1, one finds 
( see sketch) : 



and 


Further, 


v x = - cos 6 2 + cos 0 3 

2 jtr 2 2jtr 3 


w x = - - sin 62 - sin 63 

2jtr 4 2itr 2 2jtr 3 

from the geometry of the system, it is seen that 


(17) 


COS 82 


sin 8 c 


J l,2 


COS 8 3 


Zl,2 


z l,2 2 + ( f - 2yi) 2> r 3 Z l,2 2 + f2 


f - 2yi 


sin 63 


•i,2 2 + ( f - 2yi) 2 ’ r £ 


z 2 + f 2 

z i, 2 + 1 


and 


r 4 = 2y x 
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The velocity components v 1 and v 1 can now be expressed in the following 
form: 


2fz n 


v, = - 


Wi = - 




,[z 1)2 2 + (f-4!y x ) a ][* 1#a 2 + f 2 ] 

_r f (f g + zx^ 2 ) 2 - 8f yi 2 (f - yi ) 

4* l yi (f 2 + 2 lf2 2 )[z 1}2 s + (f - 2y 1 ) 2 ] 


(18) 


If equations (l4) and (l8) are used to express v x and \r 1 as functions 
of y x , then equation (l6) can be integrated directly to give the path 
of vortex 1. After algebraic manipulation, equation (l6) thus becomes 


f 4 +2f 3 yi ( G -1 ) -Gf 2 yi 2 ( G+2 ) +2G 2 f y x 3 -G 2 y x 4 


zi-z 0 = 


dyi 


y 0 (^-Gfy-L+Gy-L 2 ) -(G+4) yi 2 +f (G+4) yi -f 2 ] [Gy 1 2 -Gfy 1 +f 2 ] 


(19) 


where z G is the initial value of . 

If now the numerator of the integrand is divided by the rational 
factor in the denominator, and the substitution 

is made, then the quantity (z x - z G ) can be expressed as the sum of four 
integrals 


Zjl-Zq - 


y ” 

‘12 


-G 


Ti g dri 


y - f - 

y i 2 


y G (G + 4) L y -- y(R 2 -Ti 2 )(Q 2 +n 2 ) 


- q 2 g 




r O 2 


r O 2 


t J{ R 2 -T1 2 )(Q 2 +ti 2 ) 


,4 Jl-5 


dr] 


2f 3 p 1 " 2 


q dq 


v v<t _rl i 1 ** v- 1 -' ‘i / \-» • 1 / - y _ 

y 0 2 °0 2 


G "t 4 (Q 2 +t) 2 ) */(R 2 -ti 2 )(Q 2 +ti 2 ) G J y (Q 2 +Tl 2 )y(R 2 -Tl 2 )(Q 2 +Tl 2 )_ 

( 20 ) 


f^G 


2 J- u 2 

where R = — — and Q = 


f 2 (4-G) 


4(g+4) 




4g 
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Equation (20) is now in a convenient form if G is less than 4. This 
corresponds to the range for periodic motion. On the other hand, if G 
is greater than 4 (aperiodic), a new constant 


"5cT 


is introduced so that equation (20) becomes 


z i _z o - 
f 

4 p* 2 


•J G(G+4) 


ji- £ 


-G/ 




yi-| 




f •/ (R 2 -T1 )(ti 2 -Q 2 ) 


+ Q 2 G 


dr) 


y 0 -2 


f /TFT £ )(ti 2 -Q 2 ) + 


G 


dTj 


2f 


,3 rfi 2 


r] drj 


4 (tiM 2 ) ^-^(^-q 2 ) g , .f (ii 2 4 2 )ySv)(7-Q 2 ) 

^ o 5 


r o 2 


( 21 ) 


The third possibility is that G is equal to 4, in which case 
Q 2 = Q 2 = 0. Equations (20) and (21) then both reduce to 


f 

, y i"S 


z i" z o = “ 


n an + f 4 P 1 2 dT ) 


y - f - 

yi 2 


dq 


j!>J f yR 2 -q 2 16/2 J f qVR 2 --q 2 qM f n 2 >/W!f 

v -- v ‘ y Q ““ * 1 


r O 2 


y 0 -i 


( 22 ) 


The three cases corresponding to equations (20), (21), and (22) will 
be treated separately since the integrations will be somewhat different 
for each. 


For cases where the vortex motion is periodic, that is, when G is 
less than 4, the determination of the vortex paths requires the evalua- 
tion of the integrals of equation (20). It is seen that the fourth 
integral of this equation is elementary while the first three are 
elliptic. Thus, while the fourth integral can be evaluated simply by 


noting that q dq = ~ d(q 2 ), the other integrals require the use of the 
elliptic transformations 2 


2 

These and subsequent transformations were obtained from a comprehensive 
table of elliptic integrals prepared by Paul Byrd, Ames Laboratory, 
MCA. 
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sn 2 Ui = sin 2 9* 



sn 2 u Q 


sin 2 



k 2 = 


R‘ 


R 2 + Q 2 


Gf 

16 


(23) 


With these transformations the integrals can he evaluated and the results 
expressed in terms of elliptic integrals of the first and second kinds. 
The final solution of equation (20) , after collection of terms and eval- 
uation of the constants, is given in nondimensional form hy the expres- 
sion 


*1 _ 


4G 

16 -g' 


E(k,9 0 )-E(k,9i) 


siii cpi cos Ti sin To cos To 


8(4-G) \Vl-k 2 sin 2 Ti \/l-k* 


sin 2 T*> 


i </g(g+4) 



sin <T>i 
-k 2 sin 2 9i 


sin 9 0 


J~T- 


k^sin 9 r 


yo 

f 


(24) 


G 

where k = - and it is recalled that z Q is equal to y Q due to the 

45° configuration of the initial locations of the vortices. (Note that 
9o lies in the second quadrant and that 9i increases positively from 
9 0 .) Thus, for the case where the motion is periodic, the motion in 
the y,z plane is completely described by equations (7), (lO), (l4) 
and (24). 


If the motion in the y,z plane is not periodic, that is, if G 
is greater than 4, then the motion of vortex 1 is described by equa- 
tion (21). Here again the first three integrals of this equation are 
elliptic, but now the required transformations are 


sn 2 u x 


sn 2 u Q 


sin 2 cp x 


sin 2 cp 0 



e R 2 -Q 2 16 

R 2 " G 2 


( 25 ) 
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It should he noted 
less than Q 2 , all 
imaginary if 


that for this case (aperiodic motion) when 
four integrals become imaginary. That is. 


yi _ 1 < Q 

f 2 f 


Z 1 


is 

is 


( 26 ) 


This condition defines the asymptotes of the vortex paths for the 
aperiodic motion • This can he seen from the relative paths shown in 
figure 1. 


If the integrals of equation (21) are now evaluated for this case, 
the resulting expression for the aperiodic path of vortex 1 is 



E(k,q>o)- E(k,9i) 


+ F(k,<Po) - F(k,cp i ) + 


2(G-4) 


^tan cp 0 •Ji- k 2 sin 2 cp 0 - tan <Pi 's/l-k 2 sin 2 ^ + 


1 /G+b f \ y Q 

2 / ~q~ l tan 9o - tan 9i J + — 


(27) 


where k = -. The motion is now completely defined by equations (7). 

G 

(10), (14), and (27). 

For the special case when G is equal to k, all the integrals are 
elementary and can readily he evaluated hy making use of the trigono- 
metric substitutions 


f -p 

y. -- y -=■ 

_ 1 2 ^ 2 
cos <Pi — ; cos cpQ = ^ 

R R 


( 28 ) 


The final solution of equation (22) is then 

z i 1 / . . \ 1 / sin Ti sin <p G 

7 - - j - sln %)* 


2 \COS^ cos cp 0 


1 cos 9o (l+ sin ^ ) 1 f sin cp x 

2 cos cp x (l+ sin <P Q ) «/2* \ cos 9i 


sin 9o 
cos cpo 


\ 1 1 


•I T 


(29) 
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Actually, it can be shown that this result is obtainable from either 
equation (24) or (27) by taking the limit as G approaches 4. It is 
interesting to note that G can assume only positive values and must 
lie between 2 and infinity as can be seen from its definition given 
with equation (l4) . Thus, the vortex paths in the y,z plane have been 
obtained for the entire range of the parameter G. Figure 2 shows the 
paths of vortices 1 and 2 for several values of G. The dashed lines 
connect corresponding positions of the two vortices at successive time 
intervals . 


Complete motions in three dimensions .- In order to complete the 
three-dimensional picture of the four vortex filaments trailing behind 
the cruciform wings, a relationship must be obtained between the points 
of the paths in the y,z plane and the distance the wing has traveled 
in the stream direction. This can be done by the relationship 


d 


Ut = U /dt = U 


r 

dy, /dt 



(30) 


Equations (l4) and (l8) can again be used to express vi as a function 
of y x so that equation (30) can be integrated. The integral thus 
becomes 


d = 


8 itGfu P 1 j fg -Qy 1 (f-y 1 ) 

^ y 0 J Gyi(f-yi) -(2y x -f ) 2 


2 / >2 

yj (f-yi) 

[ f 2 -Gy 1 (f-y 1 )]^ 


dyx 


Again using the substitution 


(31) 


Y x = T| + 


_f 

2 


equation ( 31 ) can be simplified into three somewhat different forms, 
one for each of the three regimes of the parameter G. 


For the case when G is less than 4 (periodic motion), equa- 
tion ( 31 ) can be written as 


d = - 


8nfU 




T - 2 


dr, 


r J G(G + 4) y o -| ( qV ) J(-R 2 -r\ 2 )(<f+rf ) 


(32) 


where 


R == 


f 2 G 


4(g + 4) 


and Q = 


f 2 (4 


G) 


4g 
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For the aperiodic case 

8 jrfU 

d = 

iVg(g+4) 


(G greater than 4), equation (31) 

r yi 'i ( T ~ ^ 

{ _£ (ti £ -Q 2 ) V(E 2 i S )(Ti 2 -Q 2 ) 

2 


becomes 

(33) 


where 

Q 2 = - Q 2 = 

* * 4g 


while for the special case when G is equal to 4, equations (32) 
and (33) both reduce to 


d 



since here 



= 0 . 


(34) 


The three expressions for the distance behind the wing can now be 
solved by use of the proper elliptic transformations. For G less 
than 4, the result is 


d TtUf 


128 


g(i6-g 2 ) L 


E(k,<Po)-E(k ,<P X ) 


8 

4- — 

G L 


F(k,q> 0 )-F(k,? x ) 


sin cp Q °os cp 0 _ sin cpi cos Ti 


where 


a - _RL = G 2 


k^ = 


4_G V J 1-k 2 sin 2 <P Q J 1-k 2 sin 2 9i 

-- kil 


(35) 


r 2 +q 2 16 


, sin^ <Pi = 


R‘ 


, and sir. cp 0 = 


R 2 -W| 

R 2 


For G greater than 4, the resulting expression is 
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d JtUf J 32 


f r Ig 2 -i6 


E(k,qP 0 )-E(k,9 1 ) 


G 

where 


^tan cp Q J 1-k 2 sin 2 cp Q - tan cp x J 1-k 2 sin 2 4^ ^ j- ( 36 ) 

/-(vf ) 2 


R g -Q g 16 


k g = - — fr = sin* 9i = 


R G 
For G eq.ua! to 4, 


R 2 -Q 2 


2y , and sin 2 <p 0 = 


R g -Q 2 


d 

f 


itUf 

r 


sln ’’’s _ ° ln _h.. + ]0g (n-sm 1 Pi)(-cob %X + 1 (am^-sin 


_ cos 2 


COS^ 


(l+sin 9 0 )(-cos cp x ) 2 


(37) 


where 


cos 


q>i = 


yi- « 


R 2 


and cos 44 



With the aid of equations (35) , (36), and (37), the positions of 
the four rolled-up vortices are completely described 3 for all distances 
behind the cruciform wings. Thus, for G less than 4, equations (7), 
(10), (l4), (24), and (35) completely determine the motion. The other 
cases are given by the corresponding equations. 


APPLICATION TO CRUCIFORM TRIANGULAR WING AT 45° BANK 

Determination of Initial Vortex Positions 
and Vortex Strengths 

The equations thus far developed give the motions of the rolled-up 
vortices in terms of their initial positions. In order to relate these 
motions to physical cases involving a given cruciform wing at a 
specified angle of attack, the initial positions of the four vortices 
must be determined as a function of the lift coefficient. This amounts 

3 

The evaluation of the circulation T will be discussed in the next 
section. 
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to the determination of the parameter G in terms of the physical 
quantities such as lift coefficient and aspect ratio. Therefore, this 
section will be concerned with the initial locations of the four vortices 
for a given cruciform wing (banked 45°) at a given flight condition. 

Since the vortices trailing behind the cruciform wing actually leave 
the trailing edges as two flat vortex sheets, the initial locations of 
the four rolled-up vortices are somewhat fictitious. However, if each 
wing of the cruciform is considered as a plane wing in sideslip, the 
locations of the rolled-up vortices from each wing can be calculated 
from the span loading curve. Such a calculation, then, can be considered 
as yielding the positions of the four vortices of the cruciform before 
any interaction has taken place between the vortices of the two wings, 
namely, immediately behind the trailing edges. Hence, in order to calcu- 
late the initial positions of the four vortices, the span loadings of the 
separate wings must be determined. 

In reference 9 , it was shown that the load distribution on each 
wing of a pointed low-aspect-ratio cruciform wing with no body is given 
by the expression 



where Ap is the difference in pressure between corresponding points 
on the upper and lower surfaces of the individual wing. 

Equation ( 38 ) is not valid for angles of sideslip greater than — . 

dx 

The span loading is obtained by integrating this load distribution in 
the chordwise direction; that is 



so 


vis/dx 


ds/dx 



+*3 




dx 


(39) 


IS 
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Equation ( 39 ) can be separated into two integrals, and since for trian- 
gular wings ds/dx is a constant, the resulting span loading is 



In the case of the unyawed plane wing, it is clear from symmetry 
considerations that, if the vortex sheet is to roll up into two 
vortices, all the vorticity on one side of the plane of symmetry must 
eventually become rolled up into one vortex. This leads to the conclu- 
sion that each rolled-up vortex must be of strength equal to the maximum 
circulation which in this case is the circulation in the plane of sym- 
metry. A corollary of this conclusion is that each rolled-up vortex 
contains only vortices of the same sense. This seems to be in accord- 
ance with existing knowledge of the behavior of vortex sheets as well as 
with existing theories of the cumulation of vorticity as applied to 
turbulence. For the wing in sideslip, whereas there is no symmetry, the 
wing can again be considered as composed of two segments, each having 
vorticity in only one direction and hence each producing one rolled-up 
vortex. The dividing line between two such segments is at the spanwise 
station of maximum circulation so that such a division leads to the 
conclusion that each rolled-up vortex must have a strength equal to the 
maximum circulation. Thus, the vortex strength is given by 

^ = 7 ; ( c Z c )max (^ 2 ) 

Further, the distance between the vortices is given by the fact that the 
lift impulse must always be that of the wing itself. Thus, 


L T = C L t S 
PU r ( c Z c )max 


(^ 3 ) 


From equation (43) it can be seen that the distance 2s 1 can be 
represented on the span loading curve as shown in the sketch, where the 
area under the rectangle is equal to the area under the curve. 
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c z c 



As seen from the sketch, the actual location of each of the vortices 
still requires the determination of the distances a* and a 2 » It has 
already been argued that all the vorticity on one side of the maximum 
circulation must eventually be contained in one rolled-up vortex. Thus 
the two areas can be equated on either side of the maximum; that is. 


ym 


( a i + y m )( c z c ) max “ J c i c d y 

_s o 

+ s o 

( a 2 ~ ^m) ( c Z c ) m ax = / C Z C d y 


y m 


(4k) 


(Note that y m is readily calculated by successive approximations.) 

To apply the above information to the cruciform wing at 45° bank, 
one merely needs to note that for this case 

a' = p = a sin 45° = ■—= 

v 2 


and that the initial position of vortex number 1 in the banked coordi- 
nate system is now given by y Q = - 7 =- It is further pointed out that 


f = yi + y 2 


a x + a 2 

~7T Z 
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Since the lift of the cruciform wing can now be expressed as 


L = /2 pUr (ai + a 2 ) = | pU 2 C L S (45) 

and the wing area S is equal to 4s 0 2 /A, it is seen that the circula- 
tion T is given by 

I = s 0 g f 2l\ (46) 

U ai + d-2 \ A / 

Sufficient information is now at hand to determine the initial 
positions of the four vortices as a function of the angle of attack for 
a particular aspect ratio or, more generally, as a function of a/A. 
However, since the above span loadings were obtained from low-aspect- 
ratio theory, for which Cl = ^ A a,, the results can be considered as a 

function of the parameter Cj^/tcA 2 . It can be shown that this will make 
the results more general in that they will now apply with the accuracy 
of linearized theory rather than of low-aspect-ratio theory. This may be 
seen from the fact that the loadings in linearized theory (see, e.g., 
reference 12) for triangular wings are simply those of low-aspect -ratio 
theory multiplied by a constant which depends upon the aspect ratio and 
the Mach number. 

The theoretical initial positions of vortex 1 have been calculated 
and are plotted against in figure 3* From these initial posi- 

tions, the values of G were calculated (equation (l4)), and are 
plotted in figure 4. The limiting angle of sideslip of equations ( 38 ) 
and (4l) has been exceeded somewhat, as indicated by the dashed portion 
of the curve in figure 3, in order to permit some interesting observa- 
tions regarding the indicated trends . The constants appearing in the 
equations for the vortex paths have now been completely determined and 
the three-dimensional paths of the four vortices can be calculated for 
any lift coefficient and aspect ratio. An illustration of such paths 
for a typical case is presented in figure 5* 

If one recalls the types of motion associated with the various 
regimes of G, figure 3 takes on added significance. The passing from 
periodic motion to aperiodic motion simply indicates the inward move- 
ment of the initial position of the upper two vortices with angle of 
attack until their velocity in the z direction is so great that once 
they pass through the lower vortices, the latter never catch up to 
complete the cycle. If this inward movement were to continue as the 
angle of attack increases (as shown in fig. 3) then at = 0.244 

the two upper vortices would have coalesced, leaving only the two lower 
vortices which would then travel parallel to the z axis as in the case 
of a plane wing. 
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Leapfrog Distance 


As an indication of the practical importance of the motion of the 
four rolled -up vortices in the calculation of downwash, a distance 
behind the wing which is characteristic of the motion will be calculated. 
This distance is taken as the distance at which the four vortices have 
become collinear, that is, the distance behind tjtie .wing at which the 
upper two vortices are just passing through the lower two. This will be 
referred to as the leapfrog distance d L measured from the wing trailing 
edges. 

The leapfrog distance can be calculated directly from equation ( 31 ) 
by taking the upper limit of integration to be the value of y x in the 
collinear configuration* This value of y x is obtained from equa- 
tion (l^) by setting z 1 2 equal to zero. Solution of this equation 
then gives the limit of integration y x as a function of G only and 
therefore of C^/itA 2 . The definite integral of equation (31) has been 
evaluated by the methods discussed previously for the entire range of G 
and the resulting leapfrog distances are plotted (in terms of chord 
lengths) against C^/itA 2 in figure 6 . For purposes of comparison, as 
well as to provide a measure of the usefulness of the assumption that 
the four vortices are fully rolled up at the trailing edge, the distance 
for the vortex sheets to roll up (as calculated by Kaden, reference l) is 
also shown in this figure. Since the distance to roll up as shown in 
figure 6 is that for a plane wing, it should be taken as an indication 
of the average distance for the four vortex sheets to roll up. Actually, 
because of the asymmetric span loading, the lower vortices will roll up 
somewhat faster, while those from the upper wing panels will roll up 
more slowly. 

The outstanding point to be noted from figure 6 is that the phenom- 
enon of leapfrogging can occur within a few chord lengths at reasonable 
lift coefficients for low-aspect-ratio triangular cruciform wings. 

The calculations are expected to represent the physical phenomenon most 
accurately when the leapfrog distance is considerably greater than the 
distance for the vortex sheets to roll up. The curves presented in 
figure 6 have been terminated at a value of C^/jrA 2 of 0 . 176 , since 
above that value (calculated from low-aspect -ratio theory) the upper 
leading edges of the cruciform wing become trailing edges and the span 
load distributions are no longer given by equation ( 4 l) . 


Rolling Up of the Vortex Sheets 


As was mentioned at the outset of the analysis, the assumption that 
the wings of the cruciform act independently of one another is not 
strictly applicable to the present problem because once the bound 
vortices of the wings become free (i.e., immediately behind the trailing 
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edge) the two vortex sheets influence one another in their subsequent 
behavior. Thus, while the foregoing analysis has treated in some detail 
the complete motions of the fully rolled-up vortices, it has not 
considered the mutual effects of the vortex sheets during the rolling— 
up process itself. Although a complete knowledge of the rolling-up 
process could be gained only by a detailed numerical analysis, some 
important questions regarding the over-all interaction effects can be 
answered by examining the nature of the vortex sheets in their initial 
cruciform configuration. Illustrated schematically in the sketch is the 

initial vortex config- 
uration, showing the 
span loading or circu- 

z lation distribution, 

i ' which is the same on 

each wing, and the 
resulting vorticity with 
directions of rotation 
indicated by the arrows . 
In view of the arguments 
presented earlier, it is 
clear from the sketch 
that all the vortices 
contained in the dis- 
tance AO will roll up 
into a single vortex, 
and that all those in 
BC will roll up into 
another single vortex 
with opposite direction 
of rotation. Further- 
more, the strength of 
the latter vortex must 
clearly be equal to the 



maximum circulation 

r 

max* 


On the other hand, it is not clear just what will become of the 
vorticity contained in OB. If the wings did, in fact, act independ- 
ently, then OB would certainly be combined with AO to form a single 
vortex of strength r mflY . However, the presence of a plane of symmetry 
at y = 0 requires that AO and OB now be considered as separate 
vortex sheets since the vorticity at 0 vanishes due to symmetry. 

In fact, closer examination of the portion of the sheets near 0 reveals 
that, due to the directions of rotation, the broken-line sheet AOB' 
must move away from the sheet A 'OB, thus severing AB and A'B' at 0. 

It is doubtful, then, that the vorticity originally contained in 
OB and OB' will ever become entrained as a part of the upper two 
vortices. It will, in fact, later be shown, from the experiments con- 
ducted for this report, that this amount of vorticity actually forms 
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a separate sheet which is. not identified 
vortices treated thus far. 


as part of any of the four 


The foregoing observations suggest a possible improvement on the 
four-vortex analysis already carried out; that is, the strengths and 
initial positions of vortices 1 and 4 could be modified to take into 
account the fact that the vorticity in OB and OB' actually is not 
on ained in vortices 1 and 4. The amount of vorticity in OB 

is^^r re f er ^ e J J.°j as r c) and resulting modification of 
is readily calculated from the span load distribution, as shown in the 
e c . e strengths of the three vortices shown are simply given by 


c l c 



and 


r c = = S 


( c Z c )max ( c z c )o 


(,A( 


P° sitions of tha vortices are again calculated by equating 
the lift impulse before and after rolling up, but now r x is consider* 
to contain only the vorticity arising from its side of the wing; that : 

a l' ( c Z c )o = f c z c dy 


— Sr 


y m 


a c (c Z c) Q + (y m -a c )(c 2 c) max = r c 2 c dy 


( a 2~y m ) ( c z c ) max - 


+ S r 


c z c dy 


m 


2k 
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Note that a 2 is unchanged by the introduction of T c . 


With the aid of the foregoing sketch, one can again express the 
cruciform-wing lift in terms of the initial vortex positions 


L =V2 pU 


^1 ( a i' + a c ) + ^2 ( a 2 


a c) 


] - 1 pU 2 C L S 


(>i9) 


and conclude that 


£2 J~2 s 0 g f Cl 

U [7( a l' + a c) + a 2 ~ a c 1 ^ A 


(50) 


where 


7 


£1 

r 2 



The original four -vortex prob- 
lem has actually now been replaced 
by the six-vortex problem illus- 
trated in the sketch. However, due 
to the relatively small strength of 
r c , as well as to the close prox- 
imity of the two opposing vortices 
T c , the influence of these addi- 
tional vortices on the motions of 
the other four is expected to be 
small. An iterative approach to 
the exact solution can be obtained 
by neglecting this influence once 
the initial position and strength 
of have been modified to allow 

for T c . Thus the problem is again 
simplified to that of the motion of 
four vortices, but now they have 


unequal strengths. (The ratio of 
the strengths to T 2 is plotted against Cp/rtA 2 in fig. 7.) 

This fact alone complicates the mathematics to the point that the 
integrations can no longer be carried out except by numerical methods. 
Therefore, the complete motions will not be determined, but the leapfrog 
distance will be recalculated in order to illustrate the effect of the 
new assumptions on the results. The leapfrog distance is readily 
checked experimentally, as will be shown in the experimental section. 
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Modified Pour-Vortex Calculation 


The equations for the relative motion of the four vortices can 
again he set "up hy the use of the general equations (3) and (4) with a 
new statement of the center-of -gravity rule in which the vortex 
strengths are no longer taken to he equal. That is. 


yc.g. 


= Snh. . Vl + , e - constant 


zr± i'i + r 2 

If it is noted that in the initial configuration 


z x } a = yi + ya; yi = y Q ; y 2 = e (1 + 7) - ry 0 


(51) 

(52) 


then equations (3) and (4) can he used as in the original analysis, 
resulting in the following expression (in dimensionless form) for the 
relative vortex paths: 



The velocity is recalculated in the same manner as before, with 

the result 
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If this expression is now used in equation (30) to determine the 
downstream distance d, the resulting (dimensionless) form of the 
integral is 


d jtUe 

e " 2f 2 X 



It is readily seen, without substituting the expression for ^x } z/ e 
as given by equation (53) > that analytical evaluation of the integral 
of equation (55) would be extremely difficult. This integration was 
carried out numerically with the definite upper limit as given by the 
collinear configuration of the four vortices (i.e., z x 2 = 0) and the 
resulting values of the leapfrog distance d^ are plotted on 
figure 6. It is noted that the resulting curve is higher than that 
calculated in the original analysis. This is due primarily to the 
smaller strengths of the upper two vortices and to the consequently 
longer time required for them to pass through the lower two because 
of reduced induced velocities of one upon the other. 


WATER -TANK EXPERIMENTS 

f 

The motions of the vortex sheets behind a triangular cruciform 
wing of aspect ratio 2 were observed experimentally by means of a 
water tank (fig. 8). The model was mounted on a vertical track which 
was driven at uniform speed into the tank, while photographs of the 
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water surface were recorded by a motion picture camera. The traces of 
the vortex sheets were made visible by applying fine aluminum powder to 
the trailing edges of the model. The model had a span of 8 inches and 
was made of 0 . 050 -inch sheet metal. 

In general, the purpose of the water-tank studies was simply to 
illustrate the motions analyzed in the earlier portions of this report. 
(A typical series of photographs enlarged from the moving pictures is 
presented 3n fig. 9-) However, it was a simple matter to obtain 
experimental values of the leapfrog distance by means of a tape which 
moved with the model track and recorded on the movie film the distance 
traveled by the wing. The results of such experiments are presented in 
figure 6 for several angles of attack and it is observed that the agree- 
ment with the modified four-vortex calculation is remarkably good. 

The formation of a separate vortex sheet arising from the center 
portion of the cruciform wing, as discussed in a previous section, can 
be seen in the photographs of figure 9- This is the phenomenon that 
has been taken into account in the modified theory. 

It is interesting to note (fig. 9) that the two lower vortices 
extend downstream in nearly the free -stream direction as indicated by 
their positions relative to the wing-tip markers. This was also true 
for the calculated paths of figure 5 . 

An Important feature of the actual flow field as distinguished from 
the simplified model used for the analysis is the persistence of the 
vortex sheets between the rolled-up vortices as seen in figure 9 - 
Because of this fact, the theoretical paths of the vortices are not 
expected to be accurate at distances behind the wing greater than d^, 
since vortices 1 and k will then begin to become entrained in the outer 
windings of the sheet that constitutes vortices 2 and 3- However, the 
magnitude of this effect will depend upon the rate of rolling -up of 
the vortex sheets and it would be expected to be most serious at the 
lower lift coefficients (or higher aspect ratios) where the motions 
being considered are not of practical concern. 


CONCLUDING REMARKS 


A detailed analysis has been made of the motions of the four fully 
rolled-up vortices trailing behind a cruciform wing which is banked 45° • 
Equations have been developed describing the paths of the four vortices 
in three dimensions, and calculations have been made of the distance 
behind the wing at which the upper two vortices leapfrog through the 
lower two. The latter calculations were confined to wings of triangular 
plan form. The simplified analysis presented in the early portions of 
the report has been modified to account for the fact that all the 
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vorticity shed from the cruciform wing is not actually contained in the 
four rolled -up vortices. 

It was found that the leapfrog distance decreased with increasing 
lift coefficient and increased with increasing aspect ratio. For angles 
of attack up to a certain critical value of C l /itA 2 , the vortex motion is 
periodic with downstream distance, while above the critical value the 
motion is aperiodic. 

From the fact that the leapfrog phenomenon can occur within two 
chord lengths of the wing trailing edges, it is clear that downwash 
calculations even at distances less than one chord length behind 
low-aspect-ratio cruciform wings must take into account the vortex 
motions considered in this report. Once the positions of the vortices 
are known, there are methods available for approximating the downwash, 
and corrections can be made for the viscous vortex cores by assuming 
that they rotate as solid bodies . 

For angles of bank other than 45°, it is doubtful, due to the lack 
of symmetry, that calculations of the type presented here could be made 
by other than numerical procedures. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Oct. 5 , 1951. 
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Figure 2.- Vortex paths in the yz plane. 


5K 


NACA TN 2605 


33 



initial position of vortex / with tiff 



NACA TN 2605 


34 
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Figure 4.- Variation of parameter G with lift 
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Figure 



5*- Typical three-dimensional paths of four vortices trailing behind cruciform wing banked 45°, 

aspect ratio 2. Cp = 1.4l. * 
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Figure 7.~ Variation of relative vortex strengths 

with lift coefficient. 
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Figure 8.- Water tank with cruciform-wing model. 
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d/c = 0.11 


d/c = 0.9^ 





d/c = 0.33 


d/c = 1.38 




d/c = 0.61 


d/c = 1.79 


(a) Stations 1 to 6. 


,NACA^ 
A- 16567 


Figure Photographs of wake at various stations behind a triangular 
cruciform wing of aspect ratio 2. Cp 2 ? 0.66. 
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d/c = 2.24 


d/c = 3.65 



d/c = 2.83 


d/c = 4.26 






d/c = 4.81 

(b) Stations 7 to 11. 
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Figure 9»- Concluded. 
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